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Abstract 

An electromagnetic pair-creation cascade seeded by an electron or a photon in an intense plane 
wave interacts in a complicated way with the external field. Many simulations neglect the vector na- 
ture of photons by including their interaction using unpolarised cross-sections. After deriving rates 
for the tree-level processes of non-linear Compton scattering and pair creation with an arbitrary 
linearly-polarised photon in a constant-crossed field, we present results of numerical simulations 
that include the photon's vector nature. The simulations of seed electrons in a rotating electric field 
of optical frequency on the one hand support the approximation of using unpolarised cross-sections 
for tree-level processes to within around 5%, whereas on the other, indicate after 1/2tt periods, 
more than a factor 3 difference in the distribution of photon polarisations and more than a 40% 
difference in the number of pairs created by these polarisations. Moreover, these figures were seen 
to increase with simulation time. The results are particularly of relevance when the external field 
is anisotropic and the polarisation of the photon is influenced by its environment. 
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I. INTRODUCTION 



There are many examples of macroscopic phenomena originating from a repeated series of 
microscopic events. One prominent example is the process of nuclear fission, where a seed 
neutron collides with a 235 U nucleus, releasing 92 Kr, 141 Ba daughter nuclei, gamma-photons 
and other high-energy neutrons that can further propagate a chain reaction Another 
is so-called "particle showers", often used by calorimeters for detection in particle physics, 
where an incident high-energy particle is brought to radiate, e.g. by passing through matter, 
and the radiation liberates other particles which in turn can radiate and further propagate 
the shower . This can occur irrespective of whether the seed particle is charged, such 
as in the common case of electrons, or whether it is neutral in the case of photons. In such 
examples, it is typically a safe assumption that in the time between radiating or freeing 



other particles, the seed particles propagate in a simple way J4]. In contrast, the cascades of 
pair-creation and Compton-scattering events initiated in intense electromagnetic fields that 
can lead to the generation of electron-positron plasmas have a much more complicated de- 
velopment. From the moment the initial particles are created, using seeds or directly from 
vacuum, their exponential growth and recycling of the external field through absorption 
and re-emission can lead to a complicated interplay between the driving external field and 
the driven plasma, with such systems predicted to occur, for example when intense elec- 
tromagnetic fields irradiate single or collections of particles {5, 6], for example in solids 
Although electron-seeded pair creation has been demonstrated experimentally js], profuse 
positron creation with lasers has thus far been mainly demonstrated via the Bethe-Heitler 
process of pair-creation by a high-energy photon in the Coulomb field of a nucleus 0, 9]. 



In order to better understand such systems, there has been an intensification of research 



efforts to simulate such plasmas 



5K7J. Due to their complexity, to model a large number of 
particles, many approximations have to be made. The purpose of this paper is to investigate 
one such approximation, namely that the polarisation of photons propagating the cascade 
can be effectively neglected, being set to the average polarisation angle. To achieve this, 
we also present a derivation of the linearly-polarised Compton-scattering cross-section in 
a constant crossed field, which we were unable to find anywhere else in the literature, al- 
though the unpolarised cross-section and pair-creation cross-section for definite polarisation 
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have been derived some time ago 1CJ, LL3J. Most recently, arbitrary photon polarisation has 
been studied in relation to tree-level Compton scattering [12] and pair-creation [3] in finite 



laser pulses (a review of strong- field QED effects can be found in 14Hl71]). One instance 



where photon polarisation in relativistic plasmas is expected to p 
magnetic field of certain astrophysical objects such as magnetars 



ay a role is in the strong 

mm 



In the current paper, we present calculations performed for a constant-crossed-field back- 
ground as this is a good approximation when the formation lengths of processes are much 
smaller than inhomogeneities in the field. More precisely, any arbitrary, time-dependent 
background can be considered constant on the QED spacetime scale when 



£ = (e 2 p^T^ u p„) / {m 2 {>cp) 2 ) <C 1, (using the definition of £ derived in 20]) where T^ v is 
the energy-momentum tensor, x is the external-field wavevector, p is the momentum of 
particle involved, e > and m are the charge and rest-energy of a positron respectively and 
we work in a system of units in which h = c = 1. In terms of laser fields, £ is often referred 
to as the "intensity" or "classical non-linearity" parameter, £ = mxo/w, Xo = Eo/E cr , E$ 
is the electric field amplitude and E cr = m 2 /e is the critical, so-called "Schwinger" field. 
Moreover, an arbitrary, constant field can then be expressed in terms of three relativistic 
invariants: 



_ e/bggf T _ e 2 F, u F^ _ e 2 F; u F^ 

where F and F* are the electromagnetic tensor and its dual. Any function of these three 
parameters W(x, J 7 , G) can be considered ~ W(x, 0, 0), when J 7 , G «C x 2 \ 1- Such functions 
then describe processes in "crossed" fields (E ■ B = E 2 — B 2 = 0, for electric and magnetic 
field E, B, equivalently F 2 = F*F = 0). At least for laser systems, since E/E cr 1, the 
second of these inequalities is easily fulfilled and as the processes in question only become 
probable when x ^ 1> the first inequality will also be fulfilled in the current study. A peda- 
gogical description of constant crossed field Compton scattering has recently been given in 



2l|. 



The paper is organised as follows. In Sec. [TTlwe present the derivation of Compton-scattering 
of a polarised photon in a constant-crossed field, discuss the result, then present in Sec. HTT1 
the rate for creation of pairs due to an arbitrarily linearly-polarised photon, which is fol- 



3 



lowed in Sec. [TV] by a study of these two processes combined - the smallest chain of events 
considered involving a real photon that can lead to e~-seeded pair creation, and finally the 
conclusions of the theoretical sections are investigated in Sec. |V] where results are presented 
from simulations of chains of lowest-order processes to compare the effect of including po- 
larisation in pair-creation cascades. 



II. DERIVATION OF LINEARLY-POLARISED COMPTON SCATTERING RATE 
IN A CONSTANT CROSSED FIELD 




Fig. 1. One example generation in the envisaged cascade, with an on-shell photon linking Compton- 
scattering (left) and pair-creation (right) processes. Double lines represent dressed wavefunctions 
that include the interaction with the external field to all orders. 



We define a plane wave to be an external field, the vector potential of which is given by 
= A^ 1 (if), which is a function solely of the phase ip = kx. The solutions to Dirac's 
equation in such a background field for a particle of momentum p are described by the 
so-called Volkov wavefunctions 221 ] : 



1 + 



2xp 
S = —px 



Ur{p) iS. 



u r (p) 



^/2p°V ' y/2p°V 
r >o ^ K V 2(xp) 



1 + 



2xp 



(2) 
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where u r (p), u r (p) are free-electron spinors, 4 = 7 My V) V is the system volume and where 
S corresponds to the action from classical electrodynamics of an electron in a plane wave 



231 ] . The limit of a constant crossed field is then achieved by choosing A^(ip) = a^ip and 



eventually letting when all dependency on x has disappeared. The amplitude for 

Compton scattering (the left-hand diagram in Fig. [1]) is given by: 



where k is the momentum of the real photon. Employing the constant-crossed-field limit, 
choosing <p = 0, one can write this as: 



S fia = e J d 4 x e i{p2+k - pi)x F( V )- F(<p) 

e I pia p 2 a 

C2 = —' 



e 2 a 2 
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2 \ xpi xp 2 J 6 \ xpi xp 2 

By Fourier-transforming F((p) and integrating over x we acquire: 



(5) 



S fi>1 = (27r) 3 e / dr 8\p 2 + k- Pl + rx) T{ 



T(r) 



dtp F{(f) e~ irip . 



To obtain the polarised rate of Compton scattering, i? 7 , we use 

V 2 fd 3 p 2 d 3 k 



tr |5'/i,'y| ; 



(6) 



(7) 



7 2T J (2tt) 3 (2tt) 3 

where tr is the trace over spin indices, the factor 1/2 is due to an average over initial electron 
spin states and T is the system duration. Using lightfront co-ordinates for momenta p + ~ = 
(p° ±p 3 )/2, p L = (p 1 ,^ 2 ) and for co-ordinates x + ~ = x° ± x 3 , x 1 - = (x 1 ,^ 2 ) and defining 
at this point the specific co-ordinate system for the calculations with x = x°(l, 0,0,1), 
ai = (0, 1, 0, 0) and a 2 = (0, 0, 1, 0), we use the following arguments to deal with the delta- 



(8) 
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(11) 
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Eq. ©: 




l-S/vyl 2 = 


{2ixfe 2 


\Sfi,7 P = 
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(27r) 6 e 2 
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(27r) 3 e 2 



Si[{T ' ~ T P-8(r - r')r(r)rV) 



(2tt) 3 L^ 7 



5(r — r') 
dr 5 4 (Ap + rx)|r(r)| 2 



-^\Ap^)8{Ap-)\Y{r^ 



where we have defined: Ap = p 2 + k—p\ and in the final line integrated over the + component 
of the delta-function to give 



Ap^ 



x 



-((P2 + kf - m 2 ) = x ( 



p 2 k 



(13) 



2p\K ' PxX 

already using what we will shortly introduce, that all momenta are on-shell, and where we 
have defined a dimensionless interaction phase length L vr/ , following standard arguments in 
e.g. Q: 



5(r — r') 



dl 
2k 
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3 i(t — r')l 



2tt ' 



(14) 



By noting that: 

/oo 
dtp (^--y^ww), (15) 
-oo 



24j and its derivative, Ai', where, for 



can be written in terms of the Airy function Ai 
example: 

j-i r a • / \ , 2tt < rc 2 2c| r - cj/3c 3 

F l = f 1 M»); h-j^r^ ^ = -3^ + 2rg ; ^l^F^' (16) 

and _F n = d 11 ^ 1 Fi/dr 71 ^ 1 for n G N>o, which occur in the definition of T(r) in Eq. OH]), 
performing the spin trace of T(r), one arrives at: 

tr|%„ 7 | 2 _ (27r) 3 e 2 5( 2 )(A^)5(A p -) t 



tr 



|r(r,)| 2 . , 2 : > r:; 
= pi£ -ri — «eae _r 2 



+ ^xA;(|F 2 | 2 + ReF 1 F 3 ) (17) 



where e' is related to the photon polarisation e via: 

e'f ■= £ » - k»— (18) 
xk 

and Eq. (TTTj) agrees with [ijj (P. 557, Eq. (36)). The redefinition of e' in Eq. (fIH|) inspired 
by 3] is also a valid a polarisation vector, obeying e' 2 = — 1 and e'k = as required, but is 
useful in removing higher powers of k x ^ y from the spin trace. Let us use the following basis 



for the two polarisation vectors transverse to the photon wavevector (e.g. as used in 25|): 



kx < 2 - fca 1|2 ^ f E\ 

Ai, 2 = t= , Oify = Sij I — , (19) 



where E is the modulus of the electric field, then A JjAt = — 5y and Aj/c = for z, j, e {1, 2} 
as required. When xfc = the rate vanishes quicker than 1/xk, so the definition in Eq. 
(Fl9|) is sound (see also |2l) for an analysis of collinear divergences in Compton scattering). 
So for a head-on collision of photon and external- field wave- vector, Ai )2 = cti, 2 . We seek the 
rate of scattering for arbitrary linear polarisation. To this end, define the polarisation to be 
a superposition of these basis vectors 

e" = Cl A1 + c 2 A£, Cl ,c 2 GC. (20) 

Since e 2 = —1, we know c| = 1 — c 2 . When one combines the expression for the rate _R 7 with 
Eq. ffT7|) and integrates the delta-functions over p 2 , just as for the unpolarised cross-section, 



the integrand is independent of k x . However, making the observation (where momenta are 
now in units of m): 



dk„ 



Xo_Xk 



dip*, 



(21) 



where tp* is the saddle-point of the Airy functions in the problem Eq. f JT5|) . and one notes 
that this is the same interaction phase length L 9tl defined in Eq. (TH1) . then the integral 
can be performed, cancelling the L V)1 factors. The final manageable integral in k y is then 
calculated using: 



CO 

oo 



dt Ai 2 (t 2 + c) = - Aii (2) 



dt Ai(t + x) 

poo I 

/ dt Ai' 2 (t 2 + c) = [3Ai'(z) + cKAh(z)] , 

J-00 An 



—00 



r dtt 2 M 2 (t 2 +c)=-j 
j — 00 ^ 



-Ai'(z) + cAii(z) 



(22) 
(23) 
(24) 



where = 2 2 / 3 , and the derivation of the final integral was based on results from 



2q . One then arrives at the rate for Compton scattering for an arbitrarily linearly-polarised 



photon emitted in a linearly-polarised constant-crossed field: 



E [0,tt[; 
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dv 



z = fi 
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2/3. Xk 



2 cos 2 6 + 1 + 
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1 + v 2 
Xk = 



Ai'(z)+Ah(z) 
e^lVl 2 



(25) 

2 Xo k-, (26) 



Xi(xi-Xk) X\ rrr 
where we have defined a polarisation angle C\ = cos<ft, C2 = sin0 (<f) = 0,7r/2 correspond 
to polarisations A 12 respectively) and a = e 2 /4n is the fine-structure constant. As a test 
of Eq. f )25|) we note that if one defines i? 7i o = [-R 7 (0) + -R 7 (7r/2)]/2 as the rate of Compton 
scattering averaged over polarisations, the so-called "unpolarised" rate, then: 



R 



7,0 



—am 
Pi Jo 



dv 



(i + vy 



1 + 



2(1 + v 2 ) 



Ai'(z) + Aix(z; 



(27) 



which, recalling that momenta are expressed in units of m, agrees with the literature e.g. 
Q (P. 559, Eq. (49)). 



A. Linearly-polarised Compton scattering rate 

In the first plot of Fig. [21 we see how -R 7 (0) depends on the quantum nonlinearity parameter 
X as well as the polarisation angle 0, where we note the characteristic maximum around 



7 



X & 1- We define the relative difference from the unpolarised rate, Ai? 7 thus: 

— 2 ^t^ (28) 

which is plotted in the second plot of Fig. Although the relative difference is largest for 
small Xii we n ote for the optimum region around \ = 1, there still persists a maximum 
relative difference of around 50% in the generation of different polarisations. 
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Fig. 2. (Color online). Ai? 7 = 2(i? 7 (0) — i? 7j o)/(i? 7 (</>) + -R 7 ,o)i for a head-on collision of electron 
and external field wave- vector (xo = 0.01). The solid black line indicates the average polarisa- 
tion produced compared with the unpolarised differential rate (dashed line), which indicates the 
polarisation of photon that is equivalent to taking the average number produced. 



III. POLARISED PAIR CREATION RATE 

Pair-creation is a cross-channel of Compton scattering, and the polarised pair-creation rate 
can be arrived at by making the following substitution in Eq. ( fTTl) : p\ — > —p^, P2 — > P3, 
k —k [23|], where the outgoing electron and positron have momenta p% and p^ respectively. 
By following similar steps to the Compton-scattering derivation, one arrives at: 

. , —am f°° du fl r 9 ,-, . , , A , . 

"M'-wrh V»("-D ^ [ 0] W " W /' <29> 

a c rn r 2/3 Xk 4u e^/lF^p^ 2 _ 

0eO,7r; z = ^ 16 ; fi = = — ; Xs = — — »-= = ^XoP 3 ■ (30) 

(Xfe - X3) Xk m z c A 



8 



This expression can be tested in an even clearer way than i? 7 (0), by comparing R e (0) and 
R e {it/2) rates with the known rates for these polarisations, which exactly reproduce the 



expressions given in e.g. (14J. In Fig. E] we again plot how the rate depends on Xk and 
polarisation angle as well as the relative difference due to this polarisation angle. We note 
that the optimum rate for pair-creation is at a typically higher value of the quantum non- 
linearity parameter than for Compton scattering, Xk ~ 10 1,25 . Also, we see from the plot of 
Ai? e , that photon polarisations which are more likely to be produced via non-linear Compton 
scattering are less likely to lead to pair-creation and vice-versa. Due to the different shapes 
of R e and i? 7 , we will further investigate whether this compensation is seen in a cascade. 
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Fig. 3. (Color online). AR e = 2(R e ((j)) — R e fl) / {R e {4>) + R e fi) for a head-on collision and xo = 0.01. 
The solid black line represents the average photon polarisation creating pairs and the dashed line 
is the photon polarisation that represents the average rate of creating pairs. 



IV. EFFECT OF PHOTON POLARISATION ON TWO-STEP FERMION-SEEDED 
PAIR CREATION 

A. Differential rate 

Electron- (positron-) seeded pair creation in an external field e~( + ) — > e~^ + e + e~ can 
proceed via a two-step process, where the intermediate photon becomes real and then decays 



to a pair, or via a one-step process where the intermediate photon remains virtual 



23| . Until 



now, the only scenario in which the virtual process has been shown to be dominant, to the 
best of our knowledge, is by tuning the external- field frequency to exploit a resonance in 
the photon propagator in the multi-photon regime < 1) [28]. As the constant crossed 
field calculation is valid in the limit of zero external-field frequency (£ — > oo), this effect 
can be ruled out. Assuming that spin effects of the incoming and outgoing fermions are 
negligible (including spin-effects originating from radiative corrections to the Volkov states 



291]). it is supposed that taking fermion-seeded pair creation in a constant crossed field to 



be given entirely by the real process, is a good approximation. Let us write this in terms 
of the expected number of pairs per unit phase required to Compton-scatter L^ 7 and to 
pair-create L Vt£ . To do this, we use the relation mTj/p - = L tp j/(>cpj). Then the differential 
rate (probability per unit interaction phase length) for the two-step process, R je (tp), can be 
formulated as: 

HyeW) = / dv '- R e (Xk{v), 0) rT , (31) 

J dv XiXk{v) 

where Xk( v ) — Xi v /{^ + v )- One can then plot the differential rate of the electron step 
for each of the polarisation possibilities, dR 7e j/dx3, where I = is the unpolarised rate, 
I = 1,2 those corresponding to polarisation A 12 in Eq. f[T9"j) and R ie = (R je (<p))^ (where 
O,/, corresponds to taking the average over 0), displayed in Fig. |U 



The dynamics for pairs created with differently-polarised photons is very similar, as can be 
seen by the differential non-linear Compton scattering rates in Fig. HJ although the maxima 
are slightly displaced and there is a slight asymmetry between dR ie /dxz and dR ie fl/ dx% as 
compared to the fixed I = 1,2 polarisations shown by lines crossing in the plots (the plots 
are identical for X3 ~ > X4 i- e - f° r electrons and positrons). It can be seen that for a higher 
incoming x (Xi)> the proportion of energy given to the created pair becomes, on average, 
lower, as the distribution in % 3 e [0, Xi[ becomes increasingly skewed towards the lower end. 
In fact, X3 stays between 0.5 < X3 < 0.75 for 1 < xi < 10 3 . This can be explained by 
noticing that Compton scattering is most probable at Xk ~ 1 and since Xk = Xz + Xi ~ 2x3, 
that the most probable value of X3 f° r Xi > 1 (hence allowing Xk > 1) is around X3 ~ 0.5. 
One could conjecture the existence of the two types of cascade mentioned in the introduction; 
a free-particle and a field- driven cascade. For a high-% incident fermion, it would seem that 
each Compton-scattering-particle-creation event reduces the ^-factor only slightly. This is 
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Fig. 4. (Color online). Plot of the differential rate for creating the next generation of pairs, given 
an incident electron with xi = 0.5 (left) and xi = 10 (right). The solid lines for definite photon 
polarisation have slightly displaced maxima and as shown in the left-hand plot, do not always 
denote the maximum and minimum rate. 



shown in Fig. [5] where the differential cross-section of R ie in \2 (the scattered fermion) is 
plotted, and the most probable ratio of X2/X1 is marked with the solid black line. Since the 
rate R ie is expressed as a probability per unit external-field phase, although large and small 
values of xi may have the same value of R iei the probability that a pair is created in a given 
duration in the lab-frame in a homogeneous field is much higher for the higher value of x as 
it traverses more external-field phase than the lower- value x particle. Therefore Fig. |5] shows 
that if a Xi ^ 5 fermion produces a pair, further acceleration from the external field will be 
required before pair-creation becomes comparatively probable again. This would represent 
a transition from the free-particle to the field-driven cascade. We note that we have focused 
simply on the two-step pair-creation process, but the two-photon Compton scattering process 



301 ] would most likely be more important to describe the fermion dynamics, as single-photon 



Compton scattering is more probable than pair-creation for all values of x considered here. 
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Fig. 5. (Color online). The differential rate of the electron step with respect to the scattered 
electron. The solid black line marks the most probable X2/X1 ratio for the scattered electron after 
emitting a real photon that decays into a pair and the dashed line marks the resulting value (%2 
becomes xi f° r t ne next generation) . The presence of a tail around X2 = Xi hints at a cascading 
process. 



B. Total rate 

By integrating under the curves in Fig. HJ we acquire an approximation to the fermion two- 
step rate, i? 7e (</>), plotted in Fig. O To deduce the overall difference that each polarisation 
makes, we plot the relative difference, A[R] with respect to R iei A[R] := (R — R ie )/R ye , 
displayed in Fig. [6j We notice that although the total rate for photons polarised in the 
I = 1,2 direction varies with Xii t ne difference between unpolarised -R 7ej o and polarised 
R*ye rates remains small at approximately 5%, whereas the maximum difference due to 
polarisation is around a 30% reduction. This represents a smoothing out of the larger 
relative differences found for the individual processes in Figs. [21 [21 In order to investigate 
the effect of polarisation when a greater variety of chains of processes occur, we turn to 
simulation. 
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Fig. 6. (Color online). Plot of the total approximation to the two-step process R ie (left) and 
the relative difference in this approximation from using a polarised intermediate photon A[R] := 
(R - Ry e ) /Rye (right). 



V. CASCADE SIMULATION 



We wish to investigate the cumulative effect of photon polarisation when Compton-scattering 
and pair-creation processes form a cascade. To this end, we employ simulation methods de- 
veloped in {5, Q], which integrate over these lowest-order rates to approximate chains of 
events. Simulation has the veritable advantage that all possible chains of real processes 
are considered, for example that several Compton-scattering steps can occur before a pair- 
creation step, which for some values of \i have the potential to expose the polarisation 
behaviour. Moreover, although we have analytically investigated the idealised background 
of a constant crossed field, in a simulation, one can employ the so-called "locally-constant 
field approximation," in which the constant crossed field rates are integrated over the phase 
of a more complicated field. In this case, \ — > x[E{(p)], where E is the electric field ampli- 
tude with a more complicated structure and (p is its phase. We take the range of validity of 
such an approximation to be the same as the validity of the constant crossed field expression, 
already discussed in the paragraph below Eq. (CQ). However, the higher the variety of chains 
of events, the more challenging it is to directly compare with theory. 
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In order to incorporate the polarised cross-sections, we rewrite the rate equations in terms 
of rates per unit energy £: 

{aix(x) + ( 2C0S2 / +1 + XkV^j Ai\x) } , (32) 

Ah(x) + ( 2COS ^ +1 - XkV^) Ai'(*)} . (33) 

where x = /x 2//3 , and the expressions resemble those of Eqs. (2,3) in (|. As an example 
scenario, we simulate the presence of 10 3 initial electrons in two cases: beginning at rest 
(Xi(t = 0) = 0) and counter-propagating with initial xi — 5 against a rotating electric field 
E(t) = (E cos x°t, E sin x°t, 0), where x° = 1 eV is the angular frequency and £ = 10 4 
(Xo ~ 0.02) (the strong-field QED effects in this field are to a good approximation equivalent 
to those of a constant crossed field background, see also (]]). For each scenario, the four dif- 
ferent cases are simulated in which: i) the parameter is randomly selected from a uniform 
distribution e [0, 7r/2[ (the more physical case) denoted without an extra subscript; ii) 
unpolarised rates are used (equivalent to setting — 7r/4), denoted by ,o; hi) = 0, denoted 
by t i; iv) = n/2, denoted by The number of fermions with p° > 20 (x > 0.4(1 — cos#), 
where 9 is the angle between p and x) and photons with Xk > 1 generated with the above 
parameters are plotted in Figs. [8] respectively. The straight part of the plots represents 
an equilibrium between momentum change due to QED processes and acceleration by the 
field. We note that the ratio of photons to fermions is of the order of unity (see Fig. [9]), 
although the number of fermions created by >c°t = 1 shows that, depending on photon po- 
larisation, the average number of generations in the cascade is between 1.8 and 3.3 (using 
N e = 2(2™ — l)iV 7 (t = 0) as the number of fermions created after n generations). 



A. Photon sector 

In a plot of the numbers of photons generated in Fig. [TJ we make the following observations: 
i) the production of photons with polarisation = is considerably more probable than 
with = 7r/2; ii) the difference in the photons generated using unpolarised and polarised 
rates is very small (this was also reflected in the frequency spectra) and iii) the time required 
for photon-production to begin was more than two orders of magnitude larger for the Xi (t = 
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Xo=0-02, Xl {t=0)=0 





'8.0 0.2 0.4 0.6 0.8 1.0 



Fig. 7. (Color online). A plot of the logarithm of the number of hard photons (\k > 1) when 10 3 
electrons initially with xi = 0, 5 respectively, interact with a rotating electric field of frequency 
x° = 1 eV, £ = 10 4 . In the inset is a log-log plot of the initial stages of the cascade. 



0) = case. For the first point, the ratio of</> = 7r/2to0 = O photons can be verified 
by calculating the ratio of the rates of photon production in each of the two cases. It 
is not immediately clear that this should be a good approximation as a few generations 
have been created, whereas the scattering rate is just for the first generation. However, it 
appears that the single Compton-scattering formula has the good agreement shown in Fig. 
ITUl For the second point, as shown in Figs. [2J El taking the unpolarised rate is equivalent to 
taking the average rate. So the observation that taking a random polarisation of photon or 
taking the unpolarised rate makes little difference is simply indicative that after thousands 
of Compton-scattering events, using the average rate for each event is a good approximation. 
Furthermore, this is supported by noting that iV 7 and iV 7i o are at the average of the positions 
of iV 7 ]1; 2 just as for the predicted average rate. The final point about the time of onset of 
Compton scattering being larger for xi(t = 0) = is also intuitive. As xi > 1 for the other 
case, Compton scattering can proceed immediately (the time step for the simulation was 
equivalent to x°t = 10~ 4 , which is exactly when pair-creation begins in the right-hand plot 
of Fig. [7]), whereas for Xi(t = 0) = 0, the seed electrons must first be accelerated before 
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! ! ! 1 1 4.5 I ! ! 1 




8.0 0.2 0.4 0.6 0.8 1.0 °8.0 0.2 0.4 0.6 0.8 1.0 



Fig. 8. (Color online). A plot of the logarithm of the number of fermions created in the same 
set-up as in Fig. [7J When production begins, there is a jump in N e due to charge conservation 
(iV e is an even number). 

they can produce photons. In a rotating electric field, one can show that the phase required 
for one of these electrons to reach %i = 1 is x°t* = a/ x°/m/xo (see }31|). This corresponds 
to t* ~ 0.07, which is comparable with t* ~ 0.05 from the simulation in Fig. [7J 



B. Fermion sector 

For the plot of the number of created fermions in Fig. El we can make the following obser- 
vations: i) although the photon-seeded pair-production rate satisfies i? 6j2 > R e ,i, we notice 
N e ,2 < N e> i] ii) the difference in the numbers of created pairs using the unpolarised and 
polarised rate is again small although they are no longer in the middle of the maximum and 
minimum curves as was the case in the photon sector iii) the time taken for pair-creation to 
ensue in the electron ensemble initially at ensemble was orders of magnitude larger than for 
Xi{t = 0) = 5. The first point can be explained by noting that any pairs created must have 
gone through a process of iV-fold Compton scattering (N > 1) followed by photon-seeded 
pair creation. Since pair-production from any low-energy photons created by Compton scat- 
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Fig. 9. (Color online). A plot of the ratio of total pairs to created hard (xk > 1) photons g ie . 
Although pairs are created much earlier in the xi(i = 0) = 5 than in the xi(i = 0) = case, after 
just 1/27T periods of the external field, the ratios already approximately equal. 



tering is exponentially suppressed R e (Xk <1)~ Xfc e ~ 8 ^ 3xfc /k° 10, lilt] , we can surmise that 
only Xk > 1 are relevant on the simulated time scales for pair-creation (in the simulations, a 
useful approximation was implemented that only Xk > 1 photons were permitted to create 
pairs). We approximate the ratio of generated pairs from photons polarised with I — 1,2: 
p e = N e)2 /N ej i using the two-step approximation in Eq. fl3Tl) . This approximation is plotted 
in Fig. [TT] and predicts the correct range of values for around the beginning of the equilib- 
rium period. For larger numbers of Compton-scattered photons, this ratio should be even 
smaller (the smaller ratio for each Compton-scattering event would likely be amplified), and 
for larger values of xPt, this is also what is observed in the numerics. For point ii), using 
again the two-step approximation, we note that in the plot of the total rate (Fig. [6]), R ie 
and -R 7ei o are very close to and even larger than for the maximum polarisation i2 7e> i, for 
Xi < 50 (typical for most fermions in the simulation). This is reflected in both plots in Fig. 
[TlT where the case xi{t — 0) =5 clearly shows this behaviour for times before equilibrium. 
The difference between using polarised and unpolarised rates for the two-step process was 
shown in Fig. [6] to be around 5%. The difference estimated from the pair-creation rate 
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Fig. 10. (Color online). A comparison of the ratio of Compton-scattering rates for the polarisations 
<p = tt/2,0 for a range of xi relevant to the simulation (left-hand plot), with the ratio of numbers 
of photons generated in the simulation (right-hand plot). The bounds of the shaded region in the 
latter are given by the left-hand plot, and the vertical dashed lines denote when the "equilibrium" 
phase of constant exponential growth in Fig. [7] is entered. 

for 1 < xi < 100 and the estimate of zero for the difference from Compton-scattering are 
compared to relative yields from simulation p. = N. — N. /N. for photons (p 7 ) and pairs 
(p e ) in Fig. [12] which, taking into account statistical fluctuations, are around the same order 
of magnitude. 



C. Discussion 



The results of numerical simulation in the current section and theoretical analysis of Sees. 
UlrllVI are in broad agreement. Moreover, the numbers of particles created when polarised 
tree-level rates were used agreed to within 5% of when unpolarised tree-level rates were 
used. This was in line with theoretical predictions from Sec. IIVI for the two-step electron- 
/positron- seeded pair-creation process when the polarisation of the intermediate photon was 
averaged over. One can surmise that the large number of seeds in the simulation aided this 
polarisation averaging effect, which would be reduced when the number of events is small. 
It was also shown that the number of photons generated with polarisation A x could be much 
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Fig. 11. (Color online). A comparison of the ratio of Compton-scattering rates for the polarisations 
<p = it/2, and a range of xi found in the simulation (left-hand plot), with the ratio of numbers of 
photons generated in the simulation (right-hand plot), where the bounds of the shaded region are 
given by the left-hand plot, and the vertical dashed lines denote when the "equilibrium" phase of 
constant exponential growth in Fig. [7] is entered. 



larger than with polarisation A 2 and that this led to a significant difference in the number 
of pairs created by these two polarisations. To put these results in context, since the photon 
polarisation is transverse to its wavevector, and since the simulation was carried out well into 
the equilibrium region where initial seed acceleration was shown to make little difference, one 
might speculate that the photon wavevectors and hence polarisations are in general isotropic. 
However, when one takes into account the fact that those photon wavevectors in the negative 
z hemisphere are more likely to be generated as Compton scattering in this direction is more 
likely, with radiation emitted in the I/7 emission cone of a relativistic fermion experiencing 



bremsstrahlung 32|, (7 = (1 — (v/c) 2 ) where v is particle velocity and c is the speed 



of light in vacuum) (see also e.g. 33[), and that these photons are also more likely to 
lead to pair-creation, then, broadly speaking, the polarisations <fi = 0, 7r/2 correspond to the 
(0, 1, 0, 0) and (0, 0, 1, 0) directions respectively. The simulations presented then demonstrate 
the principle that, in a system where an anisotropy in the photon vectors exists (such as non- 
linear Compton-scattering from highly- relativistic electrons in an external field), polarisation 
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Fig. 12. (Color online). In the plot are the relative differences between the number of photons and 
pairs produced for the polarised and unpolarised case (p. t o = N. — N.q/N.). The solid red horizontal 
line and grey-blue zones are the predictions using rate arguments for the respective curves. 



can play an important role, particularly if this polarisation can be modified by the external 
field between events in a cascade. It is a straightforward calculation to show that if a 
photon polarised in the x-y plane with an initial angle 0o to the x-axis gains a constant 
phase change S(p x>y <C 1 along each of these axes, the subsequent rotation angle (dichroism) 
of the polarisation vector becomes Sep = — [(5(p x — 5ip y ) 2 /8] sin40 o , and the induced ellipticity 
is e = [(5ip x — 5ip y )/2] sin2eV It follows that a randomly-aligned polarisation vector will 
eventually become either ordinary or extra-ordinary depending on <po and remain so (see Fig. 
[T3l) . Likewise the ellipticity will eventually become zero and the photon linearly-polarised. 
Although an idealistic model, the simulations presented in the current paper would imply a 
modification to pair-creation rates in such an environment. 



Phase changes can be induced by the polarised vacuum, where it has been shown in more 
complicated backgrounds such as e.g. focused lasers [34j], that ip 3> e is possible, but also 
in a plasma itself, in which it has been shown that vacuum polarisation effects can also be 
enhanced 35|. Although likely irrelevant for laser-based experiments as the induced ip is 



too small to lead to a significant change in N~ 



34 



, the polarisation-dependent results 
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= 7T/2 




= 

Fig. 13. For a constant phase shift in the x- and y- directions of a photon propagating in the 
z direction, the diagram represents how the 4> = ir/4 fixed-point is repulsive and the <f> = 0,7r/2 
attractive (solid arrows) for initially random photon polarisations (dashed-line arrows). The same 
behaviour follows in the other quadrant tt/2 < <p < tt. 

derived in the current paper could be of importance in astrophysical scenarios such as in 
the field around the magnetospheres of pulsars where Yn > 0.1 is possible (one example of 
this is in the soft gamma-ray repeater SGR 1900+14, |37||). More thorough calculation and 
modelling of the strong magnetic field is required before the influence of these results in this 
area can be ascertained. 



VI. SUMMARY 

We have presented a derivation for the rate of non-linear Compton scattering and photon- 
seeded pair-creation for linearly-polarised photons in a constant crossed field. Depending 
on the fixed polarisation of the photon involved, the rates for these processes were seen to 
vary between +40% and —80% that of the unpolarised rates. Moreover, those polarisations 
of photon that were more likely to be produced by non-linear Compton scattering were less 
likely to be produced by pair-creation and vice versa. To study the combined effect of this 
disparity, the two-step electron-seeded pair-creation process in a constant crossed field was 
calculated ^ — > ^ + 7, 7 — > e + e~ . Analytical results show that when the photon has a 
fixed polarisation, for incident electron cm-parameter 1 < xi < 100 the rate for two-step 
electron-/positron- seeded pair creation can be up to 30% lower than when the photon is 
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considered unpolarised. However, the results also show that when the polarisation of the 
photon is averaged over, the difference from using unpolarised rates for each part of this 
two-step process, was only around 5%. To test whether photon polarisation plays a role 
when other chains of Compton-scattering and pair-creation events are involved, for example 
in the creation of an electron-positron plasma, we used the numerical framework developed 



in 



The results of simulations in a rotating electric field of frequency 1 eV with 
intensity-parameter £ = 10 4 for the two cases of having 10 3 initial electrons with xi — 
and Xi — 5 were shown to support these conclusions for electromagnetic cascades of on 
average between two and three generations. On the one hand, the agreement to within 
5% in the number of pairs created by using polarised and unpolarised photons was found, 
supporting this approximation when simulating electron-positron plasmas in intense lasers. 
On the other, simulations also agreed with another prediction from theory, that the difference 
between the most and least prevalent photon polarisations was more than 300%, with the 
difference in the rate of pair-creation from photons with these polarisations being more than 
40%. Furthermore, simulations hinted that both these figures grow with time. These results 
are particularly relevant when the photon vectors are anisotropic and when the photons' 
polarisation can be modified by its environment. 
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